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ABSTRACT 

This  work  deals  with  the  effect  of  reduced  pressure 
in  a  controlled  turbulent  flow  on  electric  discharge 
stabilization.  A  converging -diverging  nozzle  was  designed 
and  fabricated  to  provide  a  device  for  investigating  this 
phenomenon.   Turbulence  was  generated  by  perforated  plates 
designed  with  a  geometry  that  had  proven  to  be  successful 
in  previous  controlled  turbulent  work.   Turbulence  and 
discharge  data  are  presented  along  with  photographs  of 
the  discharge  for  no  flow,  laminar  flow  and  turbulent  flow. 
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I.   INTRODUCTION 

Several  classes  of  electric  discharge  gas  laser  devices 
are  currently  under  investigation.  (Ref .  1)   A  stagnant  or 
no  flow  device  accomplishes  heat  removal  and  discharge 
stabilization  by  diffusion  to  cavity  walls.  Maximum  prac- 
tical power  capability  for  this  class  of  lasers  is  approxi- 
mately 500  watts  (W)  with  a  10  meter  (m)  glass  tube.   Another 
one  incorporates  gas  flow  for  convective  cooling  and  dis- 
charge stabilization  by  diffusion  to  cavity  walls.   These 
devices  have  a  practical  power  capability  of  5kW  with  a 
7  m  glass  tube.   A  similar  class,  using  convective  cooling, 
increases  discharge  stabilization  by  various  pre-ionization 
techniques  with  power  capabilities  greater  than  lOkW. 

Recent  work  with  turbulent  flow  for  discharge  stabiliza- 
tion adds  a  fourth  class  of  laser  devices  which  have  prom- 
ising power  capability.   Biblarz  and  Nelson  (Ref.  2)  report 
significant  discharge  power  increases  with  ambient  pressure 
turbulent  flow.   Low  pressure  (80-200  Torr)  turbulent  flow 
studied  by  Wiegand  and  Nighan  (Ref.  3)  found  increased 
discharge  stabilization  as  a  consequence  of  the  turbulence. 
Another  low  pressure  study  by  Wasserstrom,  Crispin, 
Rom  and  Shwartz  (Ref.  4)  used  a  transverse  secondary 


flow  to  create  vortical  flow.   They  also  found  that  vortex 
flow  had  a  significant  effect  on  discharge  stabilization. 

The  power  capability  of  an  electric  discharge  connec- 
tive laser  is  scalable  with  mass  flow  rate  (Ref .  5) .   The 
mass  flow  rate  can  be  increased  by  increasing  density, 
discharge  area  and  velocity.   Elevated  pressure  causes 
difficulty  with  discharge  stabilization  but  more  power  can 
be  coupled  into  the  lasing  mixture.   No-flow  breakdown 
data  at  different  discharge  pressures  for  4  cm  and  8  cm 
electrode  gaps  are  presented  by  Loeb  in  Reference  6.   Using 
these  data,  breakdown  power  for  increasing  pressure  is 
shown  in  Figure  1.   A  single  maximum  is  reached  for  a 
given  gap  length,  indicating  an  optimum  discharge  pressure. 
The  purpose  of  this  project  is  to  investigate  the  effect 
of  sub-ambient  controlled  turbulence  flow  on  discharge 
performance  in  the  range  above  450  Torr. 

Section  II  deals  with  selecting  an  economically  fabri- 
cated device  to  explore  the  desired  region  and  discusses 
the  general  philosophy  supporting  its  design.  A  complete 
description  of  the  entire  apparatus  is  presented  in  Section 
III.   The  experiment  was  divided  into  two  primary  areas, 
turbulence  tests  and  discharge  tests,  which  are  discussed 
in  Section  IV.   The  results  of  the  experiment  and  conclusions, 
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with  some  recommendations  for  improvement,  are  found  in 
Sections  V  and  VI  respectively. 

An  interesting  model  of  space  charge  convection  effects 
is  presented  in  Reference  7.   Numerous  experiments  were 
conducted  during  this  project  to  provide  data  to  help 
establish  and  support  the  proposed  model.   Air  was  used 
for  the  discharge  medium  in  all  tests  performed  during  this 
project.   A  CCL,  N  and  He  mixture  would  be  expected  to 
produce  superior  results. 
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II.   REDUCED  PRESSURE  SYSTEM  DESIGN  CRITERIA 

The  design  for  reducing  pressure  in  the  electric  dis- 
charge test  section  resulted  from  a  preliminary  study  of 
two  gasdynamic  methods.   First,  the  possibility  of  using 
a  converging  nozzle,  such  that  the  inlet  was  open  to 
atmospheric  air  and  the  exit  coupled  to  an  evacuated 
reservoir,  was  considered.  A  pressure  ratio  of  sixty  per 
cent  of  atmospheric  air  could  be  obtained  in  the  test 
section.   However,  the  flow  rate  required  to  produce 
this  pressure  reduction  combined  with  the  need  for  at 
least  a  two  minute  run  time  dictated  the  use  of  a  reservoir 
and  vacuum  pump  much  larger  than  was  locally  available. 
This  approach,  therefore,  was  subsequently  abandoned. 

An  alternative  method  which  consists  of  using  the  flow 
characteristics  of  a  converging -diverging  nozzle  was  adopted 
The  system  is  driven  by  compressed  air  at  the  inlet,  with 
the  diffuser  exit  open  to  atmospheric  air.  A  short,  con- 
stant area  section  at  the  throat  of  the  nozzle  houses  the 
electrodes  for  the  electric  discharge  test  section.   The 
turbulence  generating  plates  used  in  the  experiment  present 
approximately  fifty  per  cent  area  blockage  to  the  flow. 
To  avoid  choking  the  flow  at  the  plate  prior  to  reaching 
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desired  test  section  conditions,  a  plate  chamber  was  located 
upstream  of  the  test  section  in  the  converging  portion  of 
the  nozzle. 

The  specific  design  of  the  nozzle  was  based  on  three 
required  features.   First,  a  pressure  ratio  of  sixty  per 
cent  of  atmospheric  air  and  Mach  0.95  were  the  desired 
maximum  operating  conditions  in  the  test  section.   Test 
section  conditions  would  be  reached  when  flow  through  the 
plate  chamber,  which  housed  a  plate  with  fifty  per  cent 
area  blockage,  reaches  Mach  0.3.   The  section  linking  the 
plate  chamber  with  the  test  section  had  to  be  as  short  as 
possible,  thereby  keeping  the  turbulence  generating  plate 
within  reasonable  range  of  the  test  section.  With  this 
configuration  discharge  density  would  be  a  function  of 
velocity  but  it  was  concluded  that  the  effect  of  these 
two  parameters  on  discharge  stability  could  be  separated. 
Details  of  the  various  sections  are  presented  in  Section  III. 

It  is  known  that  vortex  stretching  occurs  in  accelerated 
flow  of  a  wind-tunnel  contraction,  thus  attenuating  the  low 
frequency,  high  energy  turbulent  spectrum  of  the  flow.  (Ref . 
8)   The  location  of  the  turbulence  generating  plate  inher- 
ently placed  the  produced  turbulent  flow  in  this  adverse 
environment.   Low  frequency,  high  energy  turbulence  is 
known  to  have  a  stabilizing  effect  on  the  discharge.  (Ref.  2) 
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It  was  felt  that  the  plate  hole  dimensions  could  be  altered 
to  compensate  for  the  attenuation  and  that  the  plate  was 
located  close  enough  to  the  test  section  to  insure  adequate 
turbulent  flow  through  the  discharge.  This  was  later  found 
to  be  incorrect.  Turbulence  decay  was  measured  and  is  dis- 
cussed in  Section  V. 

Since  breakdown  conditions  are  slightly  variable  with 
changes  in  water  vapor  partial  pressure  (Ref.  9),  a  study 
was  conducted  to  determine  a  method  for  drying  the  air 
supply.   Chemical  dryers  provide  the  most  effective  means 
for  drying  moist  air  but  commercial  products  capable  of 
accommodating  the  expected  flow  rates  are  large,  expensive 
units  thus  precluding  their  use.   To  provide  consistency 
of  measured  data,  the  effect  of  relative  humidity  had  to 
be  taken  into  account.   It  was  decided  that  breakdown 
voltage  data  would  be  corrected  to  a  standard  humidity 
of  0.6085  inches  of  mercury  vapor  pressure  by  applying  a 
correction  factor  determined  from  curves  presented  in 
Reference  9.   The  procedure  used  for  determining  the  partial 
pressure  of  water  vapor  in  the  discharge  region,  and  the 
correction  factor  curve,  are  presented  in  Appendix  A. 
Recent  work  related  to  water  vapor  effects  on  breakdown 
voltage  can  be  found  in  Reference  10. 
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III.   EXPERIMENTAL  APPARATUS 

The  experimental  apparatus  is  composed  of  three  major 
subsystems.   These  subsystems  include  an  air  flow  system, 
a  variable  high-voltage  D.C.  discharge  circuit,  and  a 
turbulence  measurement  system.  Figure  2  is  a  photograph 
of  the  laboratory  and  most  of  the  equipment  used  in  the 
experiment. 

A.   THE  FLOW  SYSTEM 

The  flow  system  consists  of  an  air  compressor,  water- 
cooler  heat  exchanger,  flow  rate  control  valves,  a  plenum 
section  and  a  converging -diverging  nozzle.  Figure  3  is  a 
schematic  of  the  entire  flow  system.  Air  supply  is  pro- 
vided by  a  Carrier  three -stage  centrifugal  compressor  with 
a  4000  ft  /min.  capability  at  a  maximum  pressure  of  two 
atmospheres.   The  flow  is  directed  through  a  water-cooler 
heat  exchanger  that  maintains  a  stabilized  flow  temperature 
at  approximately  90°F. 

Three  gate  valves  are  used  to  control  the  flow  rate 
through  the  test  section.   Two  of  these  valves  are  located 
downstream  of  the  heat  exchanger  and  are  the  primary  control 
valves  for  the  flow.   It  was  determined  that  the  flow  choked 
in  the  heat  exchange  when  the  primary  valves  were  set  for 
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maximum  flow  and  this,  therefore,  necessitated  the  installa- 
tion of  a  heat  exchanger  by-pass  controlled  by  the  third 
gate  valve.   The  by-pass  eliminates  the  problem  and  full 
flow  is  available. 

A  plexiglass  plenum  chamber  used  in  Reference  11  was 
incorporated  into  this  apparatus  but  it  catastrophically 
failed  during  early  operation  and  was  subsequently  replaced 
with  a  3/4  in.  thick,  plywood  structure.   The  plywood 
plenum  joints  were  glued  and  reinforced.   Internal  and 
external  surfaces  were  sealed  to  prevent  moisture  absorp- 
tion.  Aluminum  honeycomb  inserts  dampen  upstream  turbulence 
and  straighten  the  flow  in  the  plenum. 

The  converging-diverging  nozzle  is  a  plexiglass  assembly 
of  five  individual  sections.   The  assembled  nozzle  is  shown 
in  Figure  4.   A  converging  section  connects  the  original 
plexiglass  plenum  chamber  to  the  plate  chamber.   The  inside 
dimensions  for  the  inlet,  i.e.,  9  in.  x  9  in.  were  selected 
to  integrate  the  new  nozzle  with  the  structure  already 
present.   The  inside  dimensions  of  the  exit  match  that  of 
the  plate  chamber. 

A  plate  chamber  houses  and  supports  the  turbulence 

generating  plates.   The  chamber's  cross  sectional  area  of 

2 
20  in.   was  selected  so  that  at  maximum  operating  conditions, 

i.e.,  Mach  0.3  flow  in  the  plate  section,  the  corresponding 
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test  section  area  (based  on  isentropic  relations,  Ref.  12) 
would  be  of  sufficient  size  to  house  the  electrodes.   Re- 
moval of  a  pressure  tight  hatch  allows  easy  access  to  the 
plate  which,  therefore,  can  be  readily  changed.   The  plate 
chamber  is  schematically  illustrated  in  Figure  5. 

A  short  converging  section  links  the  plate  chamber  with 
the  test  section.   Its  length  was  determined  by  fabrication 
factors  and  the  desire  to  have  the  plate  as  close  as 
possible  to  the  test  section.   The  test  section  provides 
support  for  the  electrodes  and  includes  a  number  of  features 
The  cathode  can  be  moved  to  provide  various  electrode 
spacings,  with  a  minimum  spacing  of  2.9  cm.   Numerous 
threaded  ports  in  the  walls  of  the  test  section  provide 
access  for  diagnostic  probes  which  are  secured  in  position 
by  pressure  tight  fittings.   The  test  section  schematic  is 
presented  in  Figure  6.   A  photograph  of  the  fabricated 
section  is  shown  in  Figure  7. 

The  diffuser  section  has  an  overall  length  of  41.18 
inches  and  a  wall  divergence  angle  of  2.5  degrees.   Tufts 
were  placed  along  the  diffuser  wall  to  check  visually  for 
flow  separation. 
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B.   THE  DISCHARGE  CIRCUIT 

Three  primary  elements  comprised  the  variable  voltage, 
D.C.  discharge  circuit.   These  elements  include  a  variable 
high-voltage,  D.C.  source,  the  discharge  electrodes  and 
current /voltage  meters.   Figure  8  is  a  schematic  of  the 
circuit. 

Two  high-voltage  power  supplies  were  available  and 
were  independently  used  in  the  experiment.   A  Sorenson 
high -voltage,  D.C.  power  supply  with  a  maximum  output  of 
30  kilovolts  (kV)  and  20  milliamperes  (mA)  powered  the 
electric  discharge  in  most  of  the  tests.   An  old  Sorenson 
Beta  high -voltage,  D.C.  power  supply  with  a  rated  maximum 
output  of  60kV  and  50mA  was  also  used  but  found  capable 
of  producing  only  34kV  and  18mA.   The  power  supplies  were 
interchanged  to  cross  check  measurements  of  experimental 
discharge  power  data.   Excellent  agreement  was  observed. 

The  discharge  electrodes  consist  of  a  pin  anode  arrange- 
ment similar  to  Aunchman's  (Ref.  11)  and  a  wire  grid  cathode, 
The  anode  pins  were  made  from  0.0625  inch  steel  wire  and 
mounted  on  an  airfoil  machined  from  0.5  to  0.125  inch  brass 
strip.   Figure  9  shows  the  details  of  the  pin  electrode. 
The  wire  grid  cathode  arrangement  resulted  from  a  series 
of  preliminary  tests  utilizing  Aunchman's  equipment.   The 
airfoil  cathode  used  by  Aunchman  was  replaced  by  an  aluminum 
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Figure  9.    Pin  Anode  Schematic 
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honeycomb  cathode  to  establish  a  more  uniform  field,  and 
did  perform  slightly  better  than  the  airfoil  arrangement. 
The  honeycomb  cathode  was  not  rigid  enough  to  withstand 
the  flow  forces  and  hence  was  replaced  by  a  wire  grid 
cathode  which  performed  just  as  well.   It  was  learned  later 
that  wire  grid  cathodes  are  the  preferred  cathode  geometries 
of  axial  flow  systems.  (Ref.  13)   A  photograph  of  the  elec- 
trodes is  shown  in  Figure  10.   Mounted  on  the  pin  anode  is 
plate  III  used  in  the  tests  and  described  later. 

Current  measurements  were  taken  with  three  Triplett 
milliammeters  with  ranges  from  0.0  to  1.0,  10.0  and  50.0 
milliamperes  wired  in  series  through  a  "make  before  break" 
switch.   A  Simpson  260  V-O-M  was  also  used  for  a  cross 
check  measurement  of  current.   Originally,  discharge  voltage 
was  measured  with  a  Sensitive  Research  electrostatic  volt- 
meter but  a  short  circuit  was  detected  in  the  meter  and 
precluded  its  use.  The  power  supply  control  panel  voltmeter 
readings  are  within  one  per  cent  of  the  Sensitive  Research 
meter  and  it  was  therefore  used  to  measure  the  discharge 
voltage.   This  error  was  taken  into  account  when  an  error 
bar  for  measured  breakdown  voltage  was  determined.   High- 
voltage  leads  were  used  throughout  the  circuit.   Power 
supply  connections  were  made  through  a  highly  polished,  oil 
submerged,  brass  sphere. 
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C.   TURBULENCE  EQUIPMENT 

Turbulence  diagnostic  equipment  includes  a  Thermo 
Systems  Inc.  hot  wire  anemometer,  a  Hewlett-Packard  oscillo- 
scope, a  General  Radio  Company,  Type  1921,  real  time 
analyzer  and  a  Hewlett-Packard  7035B  X-Y  recorder.   Hot 
wire  probe  filaments  are  0.00015  inch  diameter  tungsten 
wires  spot  welded  to  the  probe  supporting  pins.   Hot  wire 
anemometer  components  include  a  model  1051-2  monitor  and 
power  supply  and  a  model  1054A-20  anemometer  module. 

The  oscilloscope  was  used  for  setting  and  adjusting 
the  stability  and  frequency  response  of  the  anemometer. 
Difficulty  was  encountered  in  matching  probe/anemometer 
impedance  following  the  procedure  outlined  in  the  operators 
manual.   To  alleviate  the  problem,  various  lengths  of  co- 
axial cables  were  kept  on  hand  for  coarse  impedance  matching, 
followed  by  fine  adjustments  of  the  anemometer.   In  general, 
the  hot-wire  system  proved  difficult  to  use.   Considerable 
skill  and  time  were  required  to  spot  weld  the  tungsten 
wires  in  place.   The  wired  probes  were  fragile  and  required 
frequent  replacement. 

The  real  time  analyzer  is  composed  of  a  model  1925, 
30  channel  multifilter  and  a  model  1926  multichannel  RMS 
detector.   Based  on  experience  gained  in  Reference  2, 
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four -second  integration  times  were  used  during  all  measure- 
ments.  The  X-Y  recorder  enables  rapid  and  accurate  data 
collection  from  the  memory  output  of  the  analyzer.   All 
equipment  is  grounded  through  a  separate  grounding  system 
in  the  laboratory. 

D.   OTHER  DIAGNOSTIC  EQUIPMENT 

A  number  of  diagnostic  instruments  used  in  the  experi- 
ment were  fabricated  locally.   Iron-cons tantan  and  copper- 
constantan  thermocouples  were  used  for  total  temperature 
measurements.   A  pitot-static  probe,  using  either  water 
or  mercury  manometers,  provided  pressure  data.   The  thermo- 
couples and  pitot-static  probe  were  carefully  calibrated 
at  a  local  facility.   Pressure  and  temperature  measurements 
of  a  free  jet  from  Mach  0  to  Mach  1.0  were  compared  with 
highly  accurate  reference  probes.   The  calibration  pro- 
cedure resulted  in  detection  of  a  faulty  probe,  which  was 
discarded,  and  permitted  a  high  level  of  confidence  in 
experimental  measurements.   Other  diagnostic  equipment  was 
also  calibrated  according  to  operating  manuals.   Calibra- 
tion results  were  used  to  estimate  an  error  bar  for  the 
various  measured  test  parameters.   Error  analysis  is 
presented  in  Appendix  B. 
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IV.   EXPERIMENTS  CONDUCTED 

The  experiments  performed  were  divided  into  two  main 
areas.   First,  turbulence  spectrum  sampling  was  required  to 
investigate  the  amount  of  turbulence  decay  encountered  with 
accelerated  flow  in  this  apparatus  and  to  make  necessary 
modifications.   Large  eddies  must  be  present  in  the  dis- 
charge region  to  have  a  stabilizing  effect  in  breakdown 
as  evidenced  by  Reference  1  and  Reference  4.   Second,  maxi- 
mum discharge  performance  was  to  be  investigated  with 
turbulent  and  non- turbulent  flow.   The  diffuser  was  removed 
to  obtain  reference  data  at  near  atmospheric  conditions. 
With  the  diffuser  on,  density  reduction  effects  could  then 
be  studied. 

A.   TURBULENCE  TESTS 

In  Reference  11,  spectrum  distributions  were  obtained 
for  various  plates.   Plate  9  had  produced  the  most  intense, 
homogeneous  turbulence  distribution.   However,  discharge 
performance  tests  at  a  velocity  of  200  ft/sec.  (61  m/sec.) 
showed  Plate  1-6  outperformed  Plate  9.   Preliminary  tests 
using  these  plates  and  Aunchman's  test  section,  showed 
that  Plate  9  was  superior  at  higher  velocities.   Plate  9 
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geometry  was  adopted  for  use  in  the  current  apparatus. 
The  various  plate  modifications  used  are  identified  as 
Plates  I,  II  and  III.   Figure  11  describes  these  plates 
in  detail. 

A  spectrum  distribution  of  the  turbulence  generated 
by  the  pin  anode  was  recorded  at  a  test  section  velocity 
of  200  ft/sec.  (61  m/sec).   Plate  I  was  used  to  investi- 
gate spectral  decay  associated  with  flow  acceleration. 
This  was  accomplished  by  placing  Plate  I  at  the  upstream 
end  of  the  plate  chamber.  (Position  A,  Fig.  3) .   Turbulence 
spectrum  was  measured  in  the  chamber  3.5  inches  (8.89  cm) 
downstream  of  the  plate.   Plate  I  was  then  relocated  at 
the  downstream  end  of  the  chamber  (Position  B,  Fig.  3)  and 
turbulence  in  the  test  section  3.5  inches  (8.89  cm)  down- 
stream of  the  plate  was  recorded.   Comparing  the  spectra 
obtained  by  this  procedure,  acceleration  effects  on  turbu- 
lence were  determined.   Spectrum  decay  was  much  greater 
than  initially  anticipated.   Plate  II  was  manufactured  and 
tested  in  the  same  manner.   With  Plate  II  installed  at 
location  B,  spectra  in  the  test  section  were  measured  with 
and  without  the  diffuser. 
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Plate  II 

Plate  III 

I 

II 

III 

MATERIAL : 

COPPER 

PHENOLIC 

ALUMINUM 
LAMINATED  ON 
PHENOLIC 

LENGTH    (IN.) 

6.0 

5.2 

3.6 

WIDTH   (IN.) 

2.6 

2.2 

1.5 

THICKNESS    (IN.) 

0.125 

0.18 

0.25 

DIAMETER    (IN.) 

LARGE  HOLES: 

0.50 

0.375 

0.375 

SMALL  HOLES: 

0.187 

0.125 

0.125 

Figure  11.    Turbulence  Generating  Plates 
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B.   DISCHARGE  PERFORMANCE  TESTS 

This  part  of  the  experiment  involved  discharge  power 
performance  under  the  influence  of  the  various  plate  and 
nozzle  arrangements  used  in  the  turbulence  tests.   Three 
internal  configurations  were  studied  with  the  diffuser 
removed.   First,  maximum  discharge  performance  with  in- 
creasing velocity  was  measured  for  a  no-plate,  "laminar" 
condition.   Again,  maximum  discharge  performance  with  in- 
creasing velocity  was  recorded  with  Plate  II  installed  in 
the  plate  chamber.   Since  power  supply  limitations  were 
reached  prior  to  obtaining  optimum  desired  test  section 
conditions,  it  was  decided  to  manufacture  Plate  III  and 
mount  it  on  the  immediate  upstream  edge  of  the  pin  anode. 
This  configuration  is  shown  in  Figure  10  and  the  plate 
position  is  identified  with  position  C  in  Figure  3.   It 
was  anticipated  that  the  flow  would  choke  at  the  plate 
prior  to  optimum  discharge  conditions,  but  it  was  felt  that 
power  supply  limitations  would  be  reached  prior  to  choking. 
Attempts  were  made  to  obtain  a  turbulence  spectrum  of  Plate 
III  in  the  discharge  region,  but  hot  wire  tungsten  fila- 
ments could  not  survive  this  intense  turbulent  flow  even 
at  low  velocities.   All  the  above  tests  were  repeated  with 
the  diffuser  on. 


35 


A  number  of  other  tests  involving  variable  voltage 
rates  and  various  gap  lengths  were  performed  to  provide 
a  data  basis  for  a  model  presented  in  Reference  7. 
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V.   DISCUSSION  OF  RESULTS 

Spectral  data  indicate  that  a  severe  low  frequency 
turbulence  decay  occurs  for  plates  positioned  at  location 
B  (Fig.  3)  in  the  accelerated  flow.  These  data  are  pre- 
sented in  Figure  12.   Comparing  the  unaccelerated  spectrum 
with  the  accelerated  spectrum  of  Plate  I  illustrates  the 
decay  encountered.   Spectra  shifted  slightly  to  the  right 
with  increasing  velocity  as  was  observed  in  Reference  11. 
Plate  III  turbulence  spectrum  could  not  be  measured,  for 
reasons  previously  stated,  but  some  partial  data  did  show 
turbulence  intensities  above  40dB.   No  changes  in  the 
spectrum  were  observed  with  diffuser  on  and  diffuser  off 
nozzle  configurations  at  a  test  section  velocity  of 
200  ft/sec.  (61  m/sec). 

Total  mass  flow  rate  through  the  nozzle  was  computed 
from  plate  chamber  diagnostics.  Mass  flow  through  the 
discharge  region  was  less  than  total  mass  flow  because  of 
plate  and  pin  anode  blockage.   The  plate  does  not  extend 
from  wall  to  wall.   The  width  of  the  plate  is  less  than 
the  total  width  of  the  test  section  and  thus  provides  a 
"boundary  layer"  slot.   The  slot  allows  higher  mass  flow 
along  the  walls  to  remove  discharge  waste  heat  and  to 
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remove  any  water  condensation.  Mass  flow  outside  the  dis- 
charge volume  was  determined  from  diagnostics  in  that 
region  of  the  test  section  and  subtracted  from  the  total 
mass  flow  giving  the  mass  flow  rate  through  the  discharge 
volume.   Figure  13  was  constructed  from  mass  flow  rate 
data. 

An  error  analysis  for  discharge  performance  curves  is 
presented  in  Appendix  B.   The  magnitude  of  the  computed 
error  bar  was  extremely  sensitive  to  pressure  readings. 
At  high  pressures,  mercury  manometers  were  used  producing 
a  larger  error  bar  than  the  readings  from  the  water  mano- 
meters. 

Maximum  discharge  power  per  mass  flow  rate  or  energy 
per  mass  obtained  for  no  plate  "laminar"  flow  and  turbulent 
flow  using  Plate  III  located  at  position  C  (Fig.  3)  are 
presented  in  Figure  14.   These  curves  show  an  initial  de- 
crease in  energy  per  mass  to  a  minimum  and  then  rapidly 
increase.   At  low  mass  flow  rates  the  total  power  remained 
relatively  constant  causing  the  initial  decrease.   As  mass 
flow  rate  increased,  convective  cooling,  turbulence  and 
sub-ambient  density  increased  the  stability  of  the  dis- 
charge resulting  in  greater  total  discharge  power.   The 
model  presented  in  Reference  7  suggests  how  velocity, 
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turbulence  and  density  influence  discharge  stability. 
The  data  presented  here  help  support  the  validity  of  the 
proposed  model. 

Sub -ambient  density  effects  on  discharge  stability 
caused  an  increase  in  energy  per  mass  as  was  anticipated. 
The  advantage  of  turbulent,  sub -ambient  flow  can  be 
realized  from  these  curves.  At  the  power  supply  current 
limited  points,  i.e.,  the  last  data  point  for  each  curve, 
more  than  twice  the  energy  per  unit  mass  was  obtained  for 
sub-ambient  turbulent  flow  at  about  half  the  mass  flow 
rate  when  compared  with  atmospheric  "laminar"  flow. 
Another  way  of  looking  at  this  is  to  superimpose  lines  of 
constant  power  on  the  curves  of  Figure  14  as  shown  in 
Figure  15.   A  given  power  level  can  be  reached  with  sub- 
ambient  turbulent  flow  at  a  much  lower  "cost"  in  mass 
flow  rate. 

Figure  14  also  shows  that  the  laminar  curves  appear 
to  be  approaching  an  upper  bound  on  energy  per  unit  mass 
while  the  turbulent  curves  are  rapidly  rising.   Data  from 
plates  I  and  II  were  identical  to  the  laminar  curves, 
demonstrating  the  severe  turbulence  attenuation  of  the 
accelerated  flow. 


42 


1.4 


1.2- 


0) 
CO 

60 

1.0 

^ 

3 

W 

§ 

0.8- 

3 

hJ 

b 

en 

en 

1 

0.6 

od 

W 

CU 

OH 

w 

3 

0.4 

PL, 

0.2- 


0.2       0.4     0.6      0.8       1.0 

MASS  FLOW  RATE  THROUGH  DISCHARGE  (kg /sec) 
Figure  15.   Constant  Power  Curves 


43 


An  average  discharge  velocity  was  computed  from  dis- 
charge mass  flow  rates.   Convective  cooling  and  turbulent 
effects  on  discharge  stability  with  increasing  average 
discharge  velocity  can  be  seen  in  Figure  16.   Increasing 
velocity  provided  convective  cooling  which  increased  dis- 
charge stability  and  hence  greater  power  density.   Higher 
power  densities  were  obtained  with  turbulent  flow  at 
velocities  greater  than  60  m/sec.   The  slope  of  the 
"laminar"  curve  appears  to  be  decreasing  above  160  m/sec. 
Peak  power  would  probably  be  reached  at  choked  conditions. 
The  extreme  ends  of  the  curves  are  power  supply  current 
limits,  not  energy  density  limits.   Figures  17  and  18 
give  breakdown  voltage  and  breakdown  current  data  for 
increasing  discharge  velocity. 

Photographs  of  the  discharge  near  breakdown  conditions 
taken  for  no -flow,  "laminar"  flow  and  turbulent  flow  are 
presented  in  Figure  19.   The  no-flow  photograph  was  a  ten- 
second  exposure.   The  general  appearance  of  the  glow 
discharge  shows  a  low  intensity  discharge  extending  be- 
tween each  pin  and  the  cathode.   The  bright  region  at  the 
bottom  of  the  photograph  resulted  from  a  breakdown  spark 
during  the  time  exposure.   In  contrast,  the  "laminar"  and 
turbulent  flow  photographs  were  one-second  exposures.   In 
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the  "laminar"  flow  photograph,  the  individual  pin  dis- 
charges were  more  diffused.  An  intense,  short,  narrow 
glow  extending  from  the  pin  anode  tips  had  developed. 
This  is  a  region  of  high  current  density.  As  shown  in 
the  turbulent  photograph,  the  discharge  had  become  even 
more  diffused.   The  bright  region  at  the  tips  of  each  pin 
was  wider  and  visibly  more  dispersed.   The  localized 
current  density  in  this  region  has  become  more  diffused 
due  to  turbulent  mixing.   Similar  observations  are 
reported  in  Reference  3.   The  outer  or  periphery  pins  of 
the  anode  were  the  first  pins  to  develop  the  intense  bright 
region  described  above.   The  center  pins  did  not  "light 
up"  as  brightly  as  the  outer  pins. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

Even  though  maximum  operating  conditions  were  not 
reached  due  to  power  supply  current  limitations,  sub- 
ambient,  controlled  turbulent  flow  had  a  significant 
effect  on  discharge  stability  in  the  working  range  of 
this  experiment.   Turbulence  generating  plates  must  be 
kept  in  the  immediate  vicinity  of  the  discharge  (location 
C,  Fig.  3)  to  avoid  spectral  decay  in  a  wind-tunnel  con- 
traction.  Power  input  is  scalable  with  mass  flow  rate. 
The  combined  effect  of  reduced  pressure  and  turbulent  flow 
at  maximum  performance  obtained  in  this  experiment  was  to 
increase  discharge  power  by  7507o  over  ambient  laminar 
flow  at  the  same  mass  flow  rate. 

The  upper  limit  of  power  input  for  this  apparatus 
would  occur  under  choke  conditions  with  Plate  III  installed 
at  location  C,  but  a  power  supply  with  sufficient  output 
would  be  needed  to  reach  that  level.   The  slope  of  the 
turbulent  flow  curves  may  or  may  not  eventually  decrease 
in  the  same  manner  as  the  "laminar"  curves,  and  this 
needs  to  be  further  investigated.   Turbulence  spectrum 
data  could  be  enhanced  by  laser  velocimeter  measurements 
of  the  turbulent  flow. 
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At  a  given  average  discharge  velocity  the  total  dis- 
charge power  can  be  obtained  from  Figure  16  for  either 
the  turbulent  or  "laminar"  case.   This  power  is  represented 
by  a  line  of  constant  power  in  Figure  15.   The  intercepts 
of  this  constant  power  line  give  mass  flow  rates  correspond- 
ing to  the  diffuser  on  or  off  configuration.   Since  dis- 
charge area  and  velocities  are  equal,  only  density  accounts 
for  the  difference  in  mass  flow  rates.   This  conclusion 
is  supported  by  the  data. 

By  extrapolating  the  probable  path  of  the  "laminar" 
curves  in  Figure  14  (see  Fig.  20),  decreasing  density  ratio 
versus  total  discharge  power  can  be  projected.   The  ratio 
of  mass  flow  rates  corresponding  to  the  intercepts  of 
constant  power  lines  gives  p  /  p  .   The  resulting  curve  is 
shown  in  Figure  21. 

There  are  a  number  of  improvements  that  can  be  incor- 
porated into  this  apparatus.   As  mentioned  earlier,  a 
larger  power  supply  would  vastly  improve  the  capability  of 
the  system.   This  would  also  require  additional  electric 
shielding  of  the  high  voltage  circuit  to  protect  the 
operators.   The  center  row  of  anode  pins  had  a  glow  inten- 
sity less  than  the  periphery  pins  for  all  flow  and  no-flow 
tests.   This  may  have  been  influenced  by  the  cathode  wire 
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Figure  21.    Power  Versus  Density  Ratio 
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grid  frame.  Optimizing  pin  density  for  turbulent  flow  may 
also  produce  interesting  results.  An  electrode  gap  length 
of  2.9  cm  was  used  for  all  discharge  data  presented.  With 
a  larger  power  supply  the  optimum  gap  for  the  corresponding 
operating  density  of  the  discharge  region  could  be  experi- 
mentally determined  and  improve  the  performance. 

A  closed,  pressure  tight  system  as  suggested  by  Reference 
11  would  alleviate  the  problem  of  relative  humidity  and  would 
provide  independent  control  of  density  and  velocity.   A 
closed  system  would  also  enable  economic  use  of  valuable 
lasing  gas  mixtures. 
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APPENDIX  A 
PARTIAL  PRESSURE  OF  WATER  VAPOR  IN  DISCHARGE 

The  relative  humidity  and  temperature  of  ambient  air 
determines  the  partial  pressure  of  water  vapor  at  the 
compressor  inlet  by  using 

w       ?„  =  OTg 

where  P  is  the  partial  pressure  of  water  vapor,  (D  is  the 
relative  humidity  and  P„  is  the  saturation  pressure  deter- 
mined from  steam  tables.   A  detailed  steam  table  can  be 
found  in  Reference  14.   From  Dalton's  Law, 
(2)      P  =  Pw  +  Pa 

where  the  subscript  a  refers  to  the  partial  pressure  of 
dry  air.   The  specific  humidity  is  given  by 

PT7  M„  PT7  LB  OF  WATER 

(2)      r*  -£  -£  m  0.621  -E  

K    J               ?a   Ma  Pa  LB  OF  AIR 

da  a 

M  and  M  are  the  molecular  weights  of  water  and  air 

respectively. 

From  Eqs.  (1)  and  (2),  and  measured  ambient  pressure, 

P  can  be  determined,  which  in  turn,  when  substituted  into 
a 

Eq.  (3)  will  give  the  specific  humidity  of  the  flow.   The 
specific  humidity  will  remain  unchanged  provided  condensation 
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does  not  occur.   If  the  flow  temperature  drops  below  the 
local  dew  point,  condensation  will  take  place  and,  therefore, 
reduce  V  . 

To  determine  the  partial  pressure  of  water  vapor  in  the 
discharge  region,  the  static  pressure  and  temperature  of 
the  test  section  were  measured.   Assuming  for  the  present 
T   remained  constant,  Eqs.  (2)  and  (3)  are  used  to  deter- 
mine Pw  in  the  test  section. 

(4)     P. 


a     (1  +  T  /0.621) 


(5) 


w      0.621 


The  steam  tables  were  used  to  determine  P   in  the  test 
section  and  substituted  into  Eq.  (1)  to  determine  the  rela- 
tive humidity.   If  the  calculated  fl)  was  greater  than  or 
equal  to  one,  V    had  decreased  or  was  about  to  and  the 
partial  pressure  of  water  vapor  in  the  discharge  region, 
Pw,  was  equal  to  P  .   If  (D  was  less  than  1,  the  partial 
pressure  of  water  vapor  as  calculated  in  Eq.  (5)  was  correct 

A  corresponding  breakdown  correction  factor  was  then 
determined  from  Figure  22.   All  data  was  corrected  to 
standard  humidity  conditions. 
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APPENDIX  B 
ERROR  ANALYSIS  (REF.  15) 

The  uncertainty  of  the  results  depended  on  accurate 
measurement  readings  of  breakdown  voltage  and  current, 
total  flow  temperature  and  static  and  stagnation  pressures, 
The  following  is  an  estimation  of  the  uncertainty  interval 
of  these  variables. 

Breakdown  voltage: 

Vb  -  (V  +  0.2)  0.99  kV 

The  0.99  comes  from  using  the  control  panel 
voltmeter.   Breakdown  voltage  was  first 
corrected  to  standard  relative  humidity. 

Breakdown  current: 

ib  =  1  i  °*2  ^  for  °  K  Z   <  10  ^ 

I,  =  I  +  1.0  mA  for  10  <  I  <  24  mA 

D       — 

Total  Temperature: 
TQ  =  T  +  2  °R 

Stagnation  Pressure: 

P  =  p  +  0.005  PSIA  for  H20  manometers 

P  =  p  +  0.02  PSIA  for  Hg  manometers 
o   r  - 
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Static  Pressure: 

Ps  "  P  ±  0.005  PSIA  for  H20  manometers 
Ps  =  P  ±  0.02  PSIA  for  Hg  manometers 

The  intermediate  results  are  functions  of  the  test  variables 
as  stated  below. 

Power : 

P  =  VI  watts 

Flow  Temperature: 

T  =  T(TQ,  pQ,  ps)  isentropic 

Velocity: 

v  =  v(TQ,  pQ,  ps)  isentropic 

Density: 

P  =P(TQ,  pQ,  ps)  isentropic 
First,  consider  the  uncertainty  in  experimental  values  of 
power. 

1)  P  =  VI  watts 

2)  P  +   aP  =  0.99  (V  +   A  V)  (I  +   A  I) 
Subtracting  eq.  1)  from  2)  gives 

AP  =  0.99  (VI  +  V   AI  +  I   a  V  +   aV   aI)-1 
a  P  =  -0.01  VI  +  0.99  (V   a  I  +  I   A  V) 

^|  =  o.99  e^  +  ir)    -°-01 
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In  a  similar  manner  the  uncertainty  interval  was 
determined  for  the  remaining  test  quantities  and  are 
summarized  below. 


Flow  Temperature: 


T«  +   AT. 


o      *o 

"0.29 


at  -  —  -1 


N 


~       =    (    PQ+     !     APrJ    \ 
\    Ps±    I    APSI    / 


Velocity:  Q  ^ 

V-J(     1-5.08  (T0+  .  T0)     [X-(M^L)°-29]J      -X 

v        ft/sec 


Density: 


AP  2.70 


Ps 
T 


±  I    Apsl  .  ribml 

+   i  at  i  J  "lf      p    Lit*  J 


Mass  Flow  Rate: 


AM  „     A   p    +  __Av  £   ribml 

ft  p  v   '  LsecJ 
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